
PRELIMINARY INVESTIGATION OF THE EFFECT 
OF SURFACE FLUCTUATIONS ON SOUND AMPLITUDE 
IN GUIDED MODE PROPAGATION 




By 

Wesley Earl Jordan 



9 




Thtb document hub been approved fa on pubtic ae- 
leabe and bate; ttb dtbtAtbut-ion ib unlimited. 



1137070 




> 



L 

















Preliminary Investigation of the Effect 
of Surface Fluctuations on Sound Amplitude 
in Guided Mode Propagation 

by 



Wesley Earl Jordan, Jr . 
Lieutenant Commander, United States Navy 
B.S., United States Naval Academy, 1963 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN ENGINEERING ACOUSTICS 



from the 



NAVAL POSTGRADUATE SCHOOL 
December 1970 






C I 







SaL POSTGRADUATE SCHOOB 

MONTEREY, CALIF. 93940 



ABSTRACT 

A simple theoretical relation between a one-dimensional sur- 
face wave and the amplitude fluctuations of an acoustic standing 
wave in a cavity with idealized boundaries is presented. An ex- 
periment was designed to investigate the applicability of the 
theory to a more complicated configuration. A 237 cm x 117 cm 
x 9 cm water-filled cavity with pressure-release walls was excited 
in a normal mode by a source placed in one corner. The resulting 
acoustic signal was monitored at various locations throughout the 
cavity with a small probe receiver. Wind-driven surface waves 
were generated by a multi-fan blower placed at one end of the 
cavity. The relationship between the resulting fluctuation of tha 
acoustic signal and the surface fluctuation above the receiver was 
invest igatod . The results are consistent with the simple theory 
at low surface-wave frequencies. Anomalous behavior was observed 
for the higher frequencies of the surface wave spectrum. Possible 
mechanisms for this inconsistency are postulated and discussed. 
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I. 



INTRODUCTION 



In A.B. WoocUs model study of shallow water propagation [l], a 
sound amplitude fluctuation, synchronized with the passage of sur- 
face waves, is reported. Wood suggested the cause was an appre- 
ciable change in water depth in the vicinity of the source and/or 
specular scatter from the surface. The reflection and scattering 
of sound from uneven surfaces’ has been given considerable attention 
in the past twenty years [2] . Most treatments in the open litera- 
ture have involved ray theoretic approaches, treating the reflect- 
ion and scattering of plane waves from irregular surfaces, and 
have assumed no interaction with a second boundary. When a second 
boundary is present, the possibility of ducted propagation arises. 

A practical mathematical tool for investigation of surface irregu- 
larities affecting this type of propagation would be normal mode 
theory. Normal mode approaches [3,4] have been rare. 

Several experimental investigations involving normal mode pro- 
pagation have been conducted. J'.A. Scrimger [5] reported corre- 
lation between surface-wave heights and sound-amplitude fluctuat- 
ions both in magnitude and in frequency content. O.S. Tonakanov [ 6 ] 
found a linear relationship between surface-wave heights and per- 
centage fluctuation in sound amplitude and phase. 

The present work describes an experimental method for the in- 
vestigation of the effect of surface variation on sound-fluctuat- 
ion in guided mode propagation. Sound-amplitude fluctuations, 
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caused by a wind-driven surface, were measured and compared to 
the surface irregularity in the vicinity of the receiver. The 
surface was treated as a Fourier superposition of various sinu- 
soidal waves, and the experimental results at each frequency were 
compared to a theoretical prediction for sound-amplitude fluctu- 
ations based on a monofrequency sinusoidal wave surface. 



10 



II. THEORY 



A normal mode solution to the wave equation for a rectangular 
enclosure with smooth boundaries (Figure 1) can be expressed as 

$ - cos k x cos k y sin k z cos COt (2.1) 

X Z j 

where CO = circular frequency of excitation 

k^ = propagation constant for the i— - direction. 

If the boundaries perpendicular to the z axis are pressure- 

release and are located at z = 0 and z = l , then k 0 f = n 0 is 

2 3 z 3 

required. In all that follows, it is assumed that the system is 



excited at frequencies allowing only n = 1 so that the eigenvalue 

3 

kg is as small as possible. 




Figure 1. Theoretical Cavity for Guided Mode 
Propagation 

If the upper boundary is disturbed, the boundaries in the z 
direction then occur at 

z = 0 

z = i z + £F(x,y,t) (2.2) 
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where GF(x,y,t) describes the change in depth of the cavity re- 
sulting from the surface fluctuation. 

If GF(x,y jt) is small, then a solution of the form 



§ = 0 ± + e0 2 



(2.3) 



is sought where 0^ is the solution for smooth boundaries s 



hown 



above and C0^ expresses a perturbation due to the altered boundary 

at z = H . 

z 

When e is sufficiently small, the following boundary conditions 



on 0^ result [ 7 ] : 



0 . 









z=f 



F(x,y ,t) 



(2.4) 



0 , 



= 0 



z = 0 



and 0^ must satisfy the wave equation 



□ 00 = 0 . 



If the surface fluctuation is a monofrequency wave propa- 
gating across the surface z = in the x direction, F(x,y,t) may 
be written 

F(x,y,t) = $, cos(yx-f2t) ( 2 *5) 

where y = surface wave propagation constant 

Cl = circular frequency of the surface wave. 
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The solution for 0^ under these conditions is 



wher e 



0 2 ' t k 3 V° S k 2 y 



co+Q 



2 



s m 


>v 


sin 


11 iK 


s in 


“a 3 


sin 


[h Jt* 
^2 z 


s in 


5 


s in 


“3^ 


sin 


V 


sin 


[i.f 

4 z 



, 

— — cos [ (kj^+yjx- (0HO)t] + 



COS [(k +y)x+(co-0)t] + 



(i^+y) 2 - k 2 2 



( 2 . 6 ) 



(2.7) 



-<v y > 2 - k / 



(— ) 
V C ' 



(k 1+ y) 2 - k/ 



M. 2 = (^) - (ic r y) 2 - * 2 2 



If Ci « CO, then the approximation 



M-i 2 = V-* = ^-(k^) 2 -^ 2 



H, 



2 . 2 



= = k -(k^y) -k 



can be made and the complete solution simplifies to 

$ ~ cos COt cos k v[cos k,x sin k z + ekJ Acos(Qt+^))] 
2 1 J z j z 



( 2 . 8 ) 



(2.9) 



wher e 



A = 



'sin [l^z 
sin 

2 Z 



2 ^ ' x2 /sin n z \ /s in H z 

„ . + 2 ■ : —7 — ) cos 2k.x 

sin H 0 £„j (sin H^f^ylsm ' 1 




sin h 2 z 



2 z- 



13 



and 



/sm n z \ / sm p. z 

sl "( k i* r > ,t MsisTpr) sini^-yjx 



Tan?!) = - 



sin 



sm [h^z 



— ■ ,- 7 - 1 cos ( k 4 -y ) x +/ — : -— 7 — ] cos ( k -V ) x 

s in I 1 1 sin ^'2 z ' v 1 ' 1 



Notice that it is the second term in the brackets of Eq. 2.9 that 
describes the fluctuation in the acoustic amplitude and that this 
term depends only on x and z. For the subsequent analysis Eq. 2.9 
will be written in the form 

$ = cos COt cos k y [cos k x sin k z + a cos(Ot +i!) ) ] ( 2 . 10 ) 

^ J- j 



where k f - n 
3 z 



a/e = (tt/2) A. 

When k^x = tt (an antinode of the standing wave in x) , 
a/e becomes 



a/e = 



S3n |1 ^z 
sin p 

1 z , 



+ 2 



(sm |i z 



sin |J. z 



sin 11 _ A I I s in p. & 

1 z / \ 2 z , 



sin p. z 

sin 'V/ 



(2.11) 



At z = ^/2, this expression further reduces to 



a/e = tt/4 



cos 

1 z 



COS ^z 



and 



§ = cos cot cos k^V [l+a cos (fit + 0 )] . 



( 2 . 12 ) 



(2.13) 




Figure 2 
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Figure 3 



16 



On the other hand, at a node where k^x = tt/ 2 and for 



x = f^/2, Eq. 2.12 becomes 



a/e 



tt/4 



1 

COS U f 

1 2 

2 



_JL_ 

cos 'll* 

2 



and $ = a cos OJt cos k^y cos(Qt + 0) . 



(2.14) 



(2.15) 



Plots of a/e as functions of surface wave frequency for the 
two cases are presented in Figures 2 and 3« These functions dis- 
play strong dependence on both acoustic and surface-wave propa- 
gation constants. 

Notice that when y = , there is a ’’resonance 1 ' effect. The 

perturbation analysis may fail in this region, because e0 may no 
longer be small compared to 0^ . If bulk absorption of the medium 
is included, a/e no longer becomes infinite, but is still extremely 
large. For y < k^ , the surface fluctuation tends to behave as a 
rather uniform change in depth of the cavity, and the amplitude 
fluctuation behaves as if the standing wave in the z direction 
were squeezed and then released. When y > k^ , the amplitude 
fluctuation appears to react to the average fluctuation over many 
surface wavelengths, and approaches the steady-state solution 0^. 

The relation between frequency and propagation constant is 

CO/k = c ( 2 . 16 ) 

for the acoustic wave and 

fyy = v (2.17) 

for the surface wave where c and V are the respective phase speeds. 
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The surface-wave speed V is dispersive and the governing 
equation [8] is 

T y 

V 2 = (g/y + ) tanh (y^) (2.l8) 

2 

where g = acceleration resulting from gravity 980 cm/sec 

p = density 

2 

= surface tension 75 dynes/cm . 

Because of this, y varies with fi/2n as shown in Figure 4 . 

In this experiment, is of magnitude 10 rad/sec and 

nj Zj 

60 ^ 90 x 10 rad/sec. Therefore, y ~ k and 60/fi ~ 90 x 10 f , • For 

this reason, the perturbation term in Eq. 2.13 appears as a low-. 

frequency modulation of the acoustic signal, and a becomes a 

measure of the modulation amplitude. The great inequality in 

acoustic and surface wave speeds means that the surface can be 

treated as effectively "frozen" with respect to the acoustic 

processes . 

The solution for £0^ in terms of GF(x,y , t) obeys super- 
position. Thus, if the surface were a spectrum of the solution 
for £0^ would be superposition of the solutions for each component 
of Q according to Eq. 2.13. 

In particular, for a narrow band of surface wavelengths, a 
localized disturbance on the surface leads to a local disturbance 
in the acoustic standing wave. Scrimger [ 5 ] observed that the 
fluctuations in the acoustic signal were strongest when generated 
surface waves were in the vicinity of the source or receiver. This 
localization of effect is assumed to apply in this experiment. 
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III. PROCEDURE AND ANALYSIS 



A. THE EXPERIMENT 

The goal of the experimental investigation was to achieve 
quided-mode propagation in a shallow water layer, disturb one 
surface with a measurable perturbation, and observe the resultant 
fluctuations in the acoustic pressure amplitude. 

In the original concept of the experiment design, one end of a 
shallow tank with pressure-release boundaries was to be lined with 
sound absorbing mater ial and a fan manifold for generating wind- 
driven "seas" was to be suspended over the other end. By using a 
highly directive source aimed into the anechoic area, a cylindri- 
cally spreading acoustic traveling wave, which excited the n = 1 
family of modes, could be transmitted. This approach would have 
yielded the following advantages: 

1. Elimination of effects resulting from reflections 
at the vertical sides of the tank. 

2. A simple traveling wave configuration corresponding to 
propagation in the lowest mode of a cy lindr ically 
diverging wave. 

3. Limitation of the ensonified surface to a region 
excluding that in the vicinity of the fans used to 
generate the surface fluctuations. 
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It was impossible within the available time to construct a 
sufficiently anechoic tank in the frequency range of interest, so 
the method had to be abandoned. The alternative was to place a 
cyl indr ically symmetrical source in a corner of the tank and ex- 
cite a normal mode of the bounded cavity. This method raised the 
immediate problem of the effects of the highly agitated surface 
directly under the fans. (See Figure 5.) (The surface under the 
fans was beaten about by the air stream generated by the fans 
directed at the water surface.) 

The theory suggests (as discussed earlier) that it may be. 
possible to assume that for small GF(x,y ,t) , the cause of the 
acoustic amplitude fluctuation observed at the receiver is the 
surface in the immediate vicinity. Therefore, acoustic noise from 
the turbulence under the fans could be treated as a ’’distant 
source” combining incoherently with the effects generated by the 
local surface. 

Excitation of the n^ »= 1 family of modes was a primary con- 
sideration; the complications coming from model interference were 
to be avoided. To ensure efficient excitation of these modes, a 
source active over a height greater than the depth of the cavity 
was chosen. This insured that, as the depth at the source changed 
with passing waves, source strength exciting the desired normal 
mode would be smoothly changing and closely related to the local 
surface. Since sea state at the source was essentially uncor- 
related with that at the receiver in the cavity, and since the 
propagation speed of the acoustic energy far exceeded the 
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Schematic Cross Section of Experiment Set .Up 



speed of the surface waves, the effect was the equivalent of 
simply changing the amplitude of 0^ . Equation 2.13 must then 
be modified to 



§ = [l+a g cos (fit + ^ s )]cos COt cos k y[l+a cos(f2t + 0)] (3*1) 

where a is the amplitude fluctuation at the receive resulting from 
the changing conditions at the source, and lb and lb are assumed in- 
coherent. For small a and a this becomes 



$ = cos COt cos k^y[l+ + a cos(J)t + </) : ) ] 



( 3 . 2 ) 



In a similar manner, the “distant acoustic sources" contributed 
by the strong surface-fluctuations in the neighborhood of the fans 
combine incoherently at the receiver with the acoustic amplitude 
fluctuations resulting from the local surface state. Inclusion 
of these “distant sources’ 1 modifies Eq. 3.2 to 



§ - cos COt cos k y[l + ^ a 2 + a ^ + a ^ cos(f2t + 1 1 ) ] 



(3.3) 



om 



where a ^ = the amplitude^ fluctuation at the receiver resulting fr 
the distant sources. 

Thus, the sound-amplitude fluctuations resulting from the local 
surface state must be obtained from the measured fluctuations and 



\ 5 

a -\ a 
a 



+ a 



and results in 



+ 




a 




( 3 . 4 ) 



( 3 . 5 ) 



Extraction of the desired values of a and then e from the ex- 
perimental data is accomplished by the following steps: 
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1. Measure the statistics of the sea surface to obtain 
e as a function of surface wave frequency (f?); 

2. Measure the acoustic amplitude fluctuations (a ) at 

a 

the various positions in the tank; 

3. Determine the values of , and a as functions of 
surface wave frequency; 

4. Compute a by Eq. 3*5 and divide it by C at each 
frequency of interest. 

B. DESCRIPTION 

The experiment was conducted in a shallow water tank whose 
dimensions are shown in Figure 6. The tank was constructed of 
3/4 inch marine plywood set on 2 inch x 6 inch supports to pre- 
vent sagging. Water depth could be maintained uniform to within 
^ 1 mm throughout the tank. 

•A 3/8 inch neoprene'*' envelope, set in the tank and held in 
place by hydrostatic pressure, served as a pressure release 
boundary. The reflection coefficient of the neoprene was measured 
by terminating an acoustic waveguide with a sample of the material 
and conducting reflectivity tests. Values of around 0.95 for the 
reflection coefficient were obtained. The boundary condition in 
the tank was confirmed by probing the sound field and verifying 
that the pressure extrapolated to zero, within experimental accu- 
racy, as the probe was brought near a boundary. 



The neoprene rubber used is foam containing nitrogen and is 
the type used in suits for SCUBA divers. It is commonly known as 
"wetsuit” mater ial . 
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Inside Dimensions 




Tank Construction 
Figure 6 
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A "wind-driven sea" was created by five equally-spaced 
squirrel-cage blowers in a manifold at one end of the tank. The 
manifold was suspended over the tank to prevent coupling of 
mechanical vibration to the tank (Figure 7 )- Each fan had a rated 
capacity of 100 ft /min. Various seas could be generated by varying 
the angle of attack of the fan output, the height of the manifold, 
orifice size, and the number of fans (Figure 8). Greater sea 
heights were obtained by placing a l6 inch wide cover over the 
water adjacent to the fan manifold in order to hold the wind down 
on the surface (Figure 5 ). 

The wave heights were measured by an immersion probe whose 
impedance varied inversely with its submerged length. This wave 
probe consisted of two submerged electrodes with a 10 kHz cur- 
rent flowing between them. The passing surface waves modulated 
the current amplitude by changing the amount of submerged elec- 
trode. This produced a current whose envelope amplitude was a 
measure of the passing wave height. After amplification, the 
information was retrieved by filtering and detecting the amplitude 
envelope of the current. 

Wave probe construction is shown in Figure 9. The positive 
electrode was made from 20-mil nikomel wire soldered to l/l6 inch 
brazing rod. The gound return was l/l6 inch brazing rod coated 
with neoprene paint over its vertical length. The last 4 cm were 
kept bare and bent at a right angle to the vertical. The elec- 
trodes were soldered to standard banana plugs and held in place 
by a connector (Pomona Electronics Connector model #l68l) . 
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Tank with Fan Manifold in Place 
Figure 7 
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Wind Generator (Fan Manifold) Operation 



Figure 8 
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The Wave Probe 
Figure 9 
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The detection and filtering network (see Figure 10) con- 
sisted of a switching diode for half-wave rectification, two 
sucessive low-pass filters with cutoff frequencies of 159-5 Hz, 
and a blocking capacitor. The 25 M»F capacitor with the 50 k Cl 
input impedance of the following amplifier formed a high-pass 
filter with 0.126 Hz cutoff. The network was tested with an 
unmodulated 10 kHz tone which ideally would generate zero DC and 
AC outputs. The AC signal which leaked through the filter amounted 
to no more than 0.05 percent of the input amplitude, and the DC 
appearing at the output was equally insignificant. 

The relation between submerged length of the wave probe and 
current flow was non-linear, the effect being more or less severe 
depending on the positioning of the probe with respect to the 
cavity boundaries. (If the electrodes were placed too near the 
bottom, the non-linearity was aggravated.) By bending the ground 
return and adjusting the probe current level , an adequately linear 
operating range was established for the wave heights expected. 

Due to the bottom effects mentioned, a static probe cali- 
bration method was used: With the probe in place, the water level 
in the tank was varied and the output voltage measured on the 
probe side of the blocking capacitor. (It was verified by use of 
an AN/URM-26 AM modulation test set that there was no measurable 
loss in amplitude across the blocking capacitor in the frequency 
range of interest.) 
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Figure 10 
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The system fox processing wave information is shown in 
Figure 11. The wave information, contained in the envelope of the 
wave probe signal, was first amplified; this information was then 
recorded on an FM tape recorder and simultaneously sent to a signal 
correlator. The output of the correlator was graphically recorded 
on the XY plotter . 

The sound source used was a BaTiO^ cylinder, 4.5 cm in 
diameter and 10.5 cm tall. The receiver was a 1/8 inch BaTiO 

3 

probe. The sound processing system is shown in Figure 12. The 
received signal was pr e-amplif ied in a low-noise 40 db amplifier, 
bandpass filtered at carrier ^2.5 kHz and then amplifier further. 
The amplitude modulation of the acoustic signal was analyzed in 
like manner with the same detector network used in the wave probe 
system. This information was then recorded in the same manner as 
the wave height data. 

The wave and acoustic amplitude fluctuation data were analyzed 
for frequency content. The system used is shown in Figure 13* The 
original tapes, recorded in the systems previously described, were 
replayed 100 times faster and recorded on tape loops with a 
standard audio recorder. The tape loops were then replayed through 
a 0.1 octave filter as shown. 

C. PROCEDURE 

The experiment was conducted in water of 8.8 cm depth. A 

strong normal mode was found at l4 . 1 kHz, a frequency well below 

cutoff of the n = 2 family of modes. The sound field was probed 
3 
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Wave Data Processing System 
Figure 11 
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The Acoustic System 
Figure 12 




The Analog Frequency Analysis System 
Figure 13 
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in order to identify the modes nodel pressures other than zero 



were found in the tank, indicating excitation of more than one 
mode or more likely imperfections in the boundaries. The cri- 
terion for a strong mode was taken to be the average standing wave 
ratio, which in this case was approximately 60 . 

The mode was identified by measuring the lengths of the 

standing waves and estimating the wave numbers. The mode with 

n = 28 , n = 9? n = 1 appeared to be excited. This mode has a 
1 2 J 

calculated frequency of 13.6 kHz (4.25 percent less than the 
observed frequency). The vertical sound profile ver if ied # that 
the n^ = 1 family of modes was being excited (Figure l4) as 
assumed in Section II. 

Once the desired mode was excited and identified, three re- 
ceiver stations were selected at strong antinodes of the hori- 
zontal wave pattern (Figure 15). The procedure now described was 
followed at each receiver position. 

The acoustic probe was set at a depth of 4.4 cm in still water, 
and the sound carrier amplitude measured. The surface was then 
disturbed with a preselected wind condition. The following ob- 
servations were then made of the envelope of the modulated sound 
car r ier : 



1. A 500 second record was made of the mean square 

amplitude of the modulation from the output of the 
correlator on the XY plotter . 
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Figure 14 
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(VRMS) 




Selected Receiver Positions 
Figure 15 
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2. A 20 minute FM tape recording was made of the 
acoustic modulation . 

3. The autocorrelation function was computed at two 
different times by the correlator and recorded 
by the XY plotter. 

The wind generator was then stopped and the surface allowed 
to calm. The carrier amplitude was again measured to verify that 
the amplitude of the standing wave had not changed. The acoustic 
data for six different sea states at each receiver position were 
obtained in this fashion. 

The area directly under the fans of the wind generator was 
intensely agitated relative to the rest of the sea surface. In 
order to isolate this contribution to the modulation of the signal, 
the tank was modified as shown in Figure l6 . The reverse flap 
lifted the wind off the surface, and the styrofoam prevented the 
propagation of surface waves. Since the styrofoam was pressure- 
release, its presence had a negligible effect on the standing wave 
pattern . 

With the probe in position at a designated receiver station 
and the wind machine operating at a specified setting, the ampli- 
tude of the acoustic modulation was measured through a 0.1 octave 
filter and average modulation amplitude as a function of frequency 
was recorded. These data were recorded at positions M B" and M C" 
only; position M A U was obstructed by the styrofoam buffer. 
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Tank Modification for Isolation 
of Acoustic Amplitude Fluctuations due to Fans 



Once all sound data had been collected, the statistics of the 



sea surface were measured. A calibration of the wave probe was 
always conducted before a data run. With the probe in place, the 
tank was filled with water to a depth of 10 cm. The water was 
then slowly siphoned out and the DC output of the wave probe was 
recorded as a function of depth. A calibration factor was de- 
termined from this curve; Probe performance did not vary with 
position in the tank except very near the walls. 

A static calibration of the wave probe is considered valid; 
however, a dynamic check was performed by simultaneously measuring 
the height of a passing wave visually and comparing it with the 
output of the probe. The passing waves were measured on a metric 
scale, printed on heavy photographic paper and taped to the side 
of the probe holder. (The thin paper ruler caused negligible dis- 
turbance in the passing surface waves.) The output of the probe 
was observed on an oscilloscope. 

Once the probe *was calibrated, it was placed at each 
receiver position and measurements were made of the amplitudes of 
the waves generated on the surface. 

The statistics of the surface at the position of the sound 
source were also required. These data could not be taken with 
the wave probe, so the behavior of the various frequency components 
of the sea were measured as they propagated down the tank and a 
projection of the behavior of each water surface at the position 
of the acoustic source was made. This estimate was checked 
visually and appeared to be valid within experimental error. 



b2 



Analog frequency analysis of the FM tape recording was made 
by a tape loop method. The 20 minute tape records were played 
back 100 times faster and recorded on 7 second tape loops with and 
audio-magnetic tape recorder. 

These loops were then played through a 0.1 octave filter with 
slow integration time. Average amplitude as a function of fre- 
quency was recorded. The tape loop output was also played through 
a true RMS voltmeter to measure overall RMS amplitude for com- 
parison with the mean square amplitude prediction of the correlator. 

D. ANALYSIS 

The RMS wave height was computed by having the correlator set 
to compute the autocorrelation function. In this mode the cor- 
relator output is 

T 

C(T) = x f V in 2(t)dt (3 ' 6) 

J o 

which equals the mean square value of the input when T approaches 
infinity. The correlator integration time is only 20 seconds, so 
a good estimate of the mean square value was obtained by taking 
the average of a 500 second graphical record of the correlator 
output (Figure 17)* The various equipment amplification factors 
were divided out and the square root taken. 

The RMS value of the envelope amplitude of the wave probe 
signal was computed in the same manner . 

The average surface wave and acoustic amplitudes as functions 
of frequency were obtained as previously discussed. As the data 
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Example of Graphical Determination of 
Estimated Value of Mean Square Amplitude 
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were recorded, a 100 Hz tone was recorded on another channel of 
the tape. The system amplification was checked by monitoring this 
test signal at various stages in the frequency analysis. This 
check also gave assurance that the data were not being distorted 
during the processing. 

In order that a and e be computed, the gains of the equipment 
through which the data was sent had to be divided out of the 
obs erva t ions . 

Refer to Figures 11 and 12 and let 



GAIN_ = the amplifier gain 
AMP 

GAIN = the FM tape recorder 

TAP b 1 

GAIN „ = the audio tape recorder. 

I AP b ^ 

The average amplitude measured through the 0.1 octave filter 
is converted to the avercige value of £ for each frequency interval 
as follows: 



WAVE PROBE OUTPUT 



FILTER OUTPUT 



GAIN 



AMP 



GAIN, 



TAPE 1 



GAIN, 



TAPE 2 



The wave probe output can be expressed in centimeters of 
wave height by dividing by the DC calibration factor determined 
before each run: 



WAVE HEIGHT 



PROBE OUTPUT 
CAL I BRAT I ON FACTOR * 



e = 



The normalized fluctuation £ then becomes 
WAVE HEIGHT _ / FILTER OUTPUT 



DEPTH 



“V GAIN TAPE r GAIN TAPE 2, 



\h 



CAL FACTORjjDEPTH 

( 3 . 7 ) 
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FILTER OUTPUT / 1 \ 

' ““amp - CAL FAcr0R • DEP ™ (““tape 1 ““tape 2 ) 

The sound amplitude fluctuation a at each frequency of 
interest is 

_ envelope AMPLITUDE AT THAT FREQUENCY 
a CARRIER AMPLITUDE 

Since the sound information was processed with the same system 
as the surface data, a was computed as follows: 



= A VERAGE AMPLITUDE 1 ( g 

“ " <“ in amp>< caer amp> (gain tape 1> ( “™tape 2> ' 

Notice that no calibration factor is required in the calculation 
of a . 



In the ratio a/e , the factor 



GAIN TAPE 1 GAIN TAPE 2 



vanishes : 



therefore, all data were treated in arbitrary units with this 
factor included. 

The carrier amplitude by which the envelope amplitude was 
divided was an average of ‘the amplitude measured before and after 
each run. If, however, the difference in these amplitudes was 
greater than 2 percent, the run was re-executed. 

The presence of high frequency sound fluctuations, unpre- 
dicted by the theory, and the persistence of these fluctuations 
throughout the tank strongly suggested the fans as a source. The 



acoustic amplitude fluctuation due to the agitation under the fans 
was obtained for positions "B" and "C" and described above. For 
any given sea state, these data were very similar at the two 



positions. Therefore, an average (fan noise) level for each sea 
state was assumed (Figure l8) for all the receiver stations. 

The other factor to be accounted for was source modulation. 

The sea state at the source was predicted as previously described.. 
The effect of the waves on the acoustic output of the source was 
determined by raising the source a given distance out of the water 
and measuring the reduction in sound amplitude at position n B n . 
With a 10 percent reduction in the radiating surface, the sound 
amplitude at position f, B n was reduced by 6 percent. This relation 
was found to be linear within the range of interest. Therefore, 
the amplitude modulation of the source was assumed to be 

a = 0.6 e , (3.9) 

S v 

(This result was checked by shading the upper part of the 
transducer the same amount with 3/8 inch neoprene foam rubber and 
comparing the results, which agreed within experimental error.) 

The amplitude fluctuation due to the local surface effect was then 
computed with Eq. 3*5: 




As a check on the validity of the data, several quantities 
were measured by two methods and the results compared. The tape 
loops, for instance, were played through the true RMS voltmeter in 
order to measure the RMS amplitude fluctuation and wave height. 

The results, adjusted for system gain, agreed to within 10’ percent 
with the RMS values computed through the correlator. 
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The frequency distribution of the modulation data was also 
determined from the autocorrelation function plotted during each 
data run. The cosine Fourier transform of the autocorrelation 
function c(7) equals the power spectral density. The transform 
of c(T) was computed by numerical integration of 100 points taken 
from the plots described above. The resulting spectra confirmed 
the distributions obtained from the tape loops. 

E. ERROR ANALYSIS 

The major cause of error- in this experiment was the non- 
stationarity of the processes. The variation in mean square 
amplitude (see Figure 17) gives a qualitative measure of the 
degree of non-s tat ionar ity involved. 

Measuring equipment with sufficiently long integration times 
to overcome this problem was not available. Therefore, long 
"effective" integration times were achieved by visual integration 
of the meter readings or .graphical averaging as previously dis- 
cussed. Study of Figure 17 suggests that an integration time of 
500 seconds is sufficient for our purposes. 

In the frequency analysis, 15 minutes of FM tape record was 

compressed onto a 7 second tape loop (see page 33)* Visually 

averaging the meter reading on the 0.1 octave filter, provided an 

integration time in excess of 500 seconds. The estimated error in 

the value assigned to acoustic amplitude fluctuation (a ) or sur- 

a 

face height fluctuation ( G ) is 10 percent. (This estimate was 
determined by comparing the values estimated by two independent 
observers . ) 
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The tape loop method was not used in obtaining the acoustic 
amplitude fluctuation due to the fans; in this case the filter was 
placed directly on the output of the filtering and detecting 
network (page 40) . Averaging two sets of observations together 
at each wind setting, minimized the error in the fan contribution 
(a^) to an estimated 15 percent, so that the more elaborate tape 
loop analysis was not deemed necessary. 

Recall that the acoustic amplitude fluctuation term a s 
resulted from the source modulation, which was extracted from 
the € curves. The estimated error in G is 10 percent. 'Adding 
a + 5 percent error for the graphical extrapolation yielding the 
surface at the source position, results in a total error in the 
source modulation term, , of +_ 15 percent. 

The RMS amplitude values computed from the graphical output 

of the correlator are considered reliable. In order to test the 

validity of the average amplitudes obtained from the tape loops, 

* 

RMS amplitudes computed in this manner were compared to those 
from the correlator. The values were consistently within 
10 percent . 

The error in a , a , and £ were assumed statistically in- 
a s f 

dependent at each receiver position. Therefore the error in the 
value of a/G at each frequency was computed using standard RMS 
error analysis. 
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IV. RESULTS AND CONCLUSIONS 



The experimental results consist of eighteen sets of data: > 

measurements were taken at three different receiver positions for 

six different winds. The six winds generated can be grouped in 

pairs: ( 1 , 2 ), (3>4) and (5? 6). Each pair of winds was generated 

with the same settings of manifold height, angle of attack, and 

orifice size of the wind machine. The even-numbered winds were 

generated with five fans; the odd-numbered ones, with three. The 

receiver positions were previously shown in Figure 15* 

The frequency distribution of the measured acoustic-amplitude 

fluctuation (a ), the fluctuation induced by the fans (a^), and the 
a x 

amplitude fluctuation due to source modulation (^ s ) are plotted 
for each wind setting at a given receiver position (see Appendix B) . 
These graphs reveal certain general features: 

1. The frequency spectrum of the measured acoustic 

fluctuation (a ) closely resembles that of the local 

surface fluctuation (e). In all cases each curve computed 

for a closely resembles the curve for the equivalent a • 

a 

2 . At each receiver position, the frequency spectra for 
both f, B !1 and n C M are similar for each wind pair. The 
stronger winds produce greater amplitudes as would 

be expected. 

3. The curves of a are similar at all receiver positions 

a 

for any given wind setting. The difference in values 

of a for positions M B" and "C" is very small, 
a 
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4. The higher frequency components of c appear to die out 

as the surface waves progress down the tank. (See 

Figure 22.) The corresponding values of a fall off 

a 

much more slowly. 

5. In computing the values of the acoustic amplitude 
fluctuation resulting from the local surface 



a 




(3-5) 



we see that a ^ becomes significant at the higher 
frequencies while a is significant at the lower 
f r equencies . 

6. There appears to be evidence of generation in a of 

a 

second harmonics of frequencies in the region corre- 
sponding to resonance. 

The computed values of a/G (using Eq. 2.13) are plotted as 
functions of surface-wave frequency for each wind condition at 
each receiver position. (Appendix C.) The results are summarized 
in Figure 19 where the curves are collected on single graphs 
according to wind condition. 

In the region below the predicted ’'resonance" the measured 
values of a/e compared favorably with the theoretical prediction. 
Very little evidence for a resonance effect exists in the results 
of any run. Above the predicted resonance the observed values of 
a/G are consistently greater than the theoretical values. The 
absence of any peaking of a/G in the vicinity of resonance cannot 
be explained. 
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Several important factors could be contributing to the high- 
frequency behavior mentioned. Recall that the term a ^ is ex- 
tremely important in determining the acoustic amplitude fluctu- 
ation a in this region. It is suspected that perhaps the entire 

observed acoustic amplitude fluctuation a above 6 Hz can be 

a 

attributed to the surface agitation under the fans. 

In the design of the fan isolation experiment, it was assumed 

that the surface agitatign directly under the fans was responsible 

for the additional amplitude fluctuation at the receiver. The 

styrofoam buffer was therefore placed under the orifice flap. 

This, however, suppressed large surface waves generated in this 

region. It is conceivable that these waves contribute to the 

tf distant sources”, but this possibility was not investigated at 

the time. Recall also that analyses of the fan-generated noise 

were not conducted with the more accurate tape- loop method 

because of time and equipment limitations. Therefore, the 

term may be larger than recorded. This hypothesis would have 

been investigated further if time had been available. 

Since accuracy in a^ becomes critical when computing the 

ratio a/e for the higher frequencies, increasingly anomolous 

behavior in a in this frequency range at distances further 
a 

removed from the fans cannot be considered as unequivocably con- 
tradicting the theory. 

The discrepancy between experiment and theory for values of 
a/e at the higher frequencies could also be attributed to 
violations of the assumptions in the theory. For instance, 
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Figure 19b 
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Figure 19c 
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Figure 20 shows that the surface fluctuations were functions of 
both x and y, rather than just of x as was assumed. Investigation 
of a more realistic expression for F(x,y,t) would require extensive 
auto- and cross-correlation on measurements of the generated sur- 
face fluctuation. 

It is also assumed in the theory that the properties of the 
surface fluctuations are uniform in the vicinity of the receiver. 
Figure 21 shows the distribution of RMS wave heights for a typical 
wind over the area of the surface in the vicinity of the receiver 
stations. This figure shows that the waves are dying out as they 
propagate down the tank. Figure 22 shows that the frequency 
spectrum for the surface waves above the frequency corresponding 
to the "resonance" is also changing. The effect of the inhomo- 
geneity of the surface could not be determined but is probably 
most severe at position "A". 

A basic assumption in the theory is that the mechanism for the 
acoustic fluctuation is the local surface state. In the investi- 
gation of the validity of this assumption, the following points 
merit study: 

1. With the normal mode excited and the water surface 

calm, the surface in the vicinity of the source was 

2 

disturbed by quadradigital agitation . An immediate 
amplitude fluctuation was observed at the distant 



Trickling four fingers in the water, making small capillary 
waves . 
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The Tank Surface Disturbed 
by a Typical Wind 

Figure 20 
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receiver (station 1! B T! ). The surface above the receiver 
was disturbed in the same manner and the acoustic signal, 
was again immediately perturbed. However, when the water 
surface was disturbed at a point removed from both source 
and receiver, an effect on the received acoustic signal 
was observed only after the surface disturbance had 
propagated into the vicinity of the source, receiver, 
and region in between. 

2. The treatment of the extreme agitation under the fans 
as a distant source (page l4) contributing to observed 
acoustic-amplitude fluctuation appears valid from the 
results of the fan isolation experiment, but this would 
appear to be inconsistent with the result just discussed. 
Recall, however, that the violent agitation of the sur- 
face under the fans is a clear violation of the assumption 
that CF(x,y,t) be small. 

The first point above appears to validate the assumption that 
the driving mechanism for the acoustic fluctuation at the receiver 
is the local surface. The fan isolation experiment- (point two) 
showed that severe surface disturbances at a remote location could 
cause amplitude fluctuations under a smooth surface. It is plau- 
sible that the distinction between these two cases results from 
the difference in magnitude of GF(x,y , t) in the two regions. 

The stringent assumptions of the simple theory were all clearly 
violated to some degree during the conduct of the experiment. It 
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is concluded, however, that the measured values of a/e below the 
frequency of predicted "resonance 11 are sufficiently close to the 
theoretical values to claim consistency between experiment and 
theory in this range. 
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APPENDIX A 



1 . 




Computer Program for a/e Calcu 
I IPLICIT REAL* 8 (A, K, P,) 


2 . 




REAL*8(DC0S, dcosh, DABS, DSQRT) 


3. 




DO 100 1 = 1 , 500 


4. 




GAM = (I-l)*.01 


5- 




Kl=2. *3. 1415926535/16. 4 


6 . 




K3=3. 1415926 5535/8 . 8 


7. 




PI~3 • 1415926535 


8 . 




AR1=2.*GAM*K1 


9. 




AR2=K3**2 


10 . 




AR3=GAM**2 


11 . 




ar4=pi/2. 


12 . 




AR5=(1 .“( (AR3+AR1)/AR2) ) 


13. 




AR6=(1.-((AR3-AR1)/AR2)) 


14. 




IF(AR5.LT.0. )GO TO 50 


15. 




AR7=DCOS(AR4*DSQRT(AR5) ) 


16 . 




GO TO 60 


17. 


50 


AR7=DCOSH(AR4*DSQRT(DABS(AR5) ) ) 


18 . 


6o 


IF(AR 6 .LT.O.) GO TO 70 


19. 




AR8=DCOS(AR4*DSQRT(AR6) ) 


20 . 


70 


AR8=DCOSH(AR4*DSQRT( DABS(AR 6 ) ) ) 


21 . 


80 


AR9=DABS(1 ./AR7-1 ./AR 8 ) 


22 . 




abe=pi*ar9/4. 


23- 




PRINT, GAM, ABE 


24. 


100 


CONTINUE 


25. 




STOP 


26 . 




END 
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APPENDIX B 



Observed Data 

The eighteen sets of data described on page 51 are plotted as 
functions of surface wave frequency (£V2 tt). 

= the measured acoustic amplitude fluctuation 
e = the normalized surface wave fluctuation in the 
vicinity of the receiver 

a ~ the acoustic amplitude fluctuation generated by the 
severe agitation under the fans 
G s = the acoustic amplitude fluctuation due to source 
modulation 

Representative error flags for these data are presented in 
Figure B1 • 
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Figure B7 
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APPENDIX C 



Measured Values of a/e 

The measured values of a/e at each receiver position for 
each wind condition are plotted as functions of surface wave 
frequency (£7/2tt). Standard RMS error flags are presented. 
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Figure C3 
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Figure C8 
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A 237 cm x 117 cm x 9 cm water-filled cavity with .pressure-release walls 
was excited in a normal mode by a source placed in one corner. The re- 
sulting acoustic signal was monitored at various locations throughout 
the cavity with a small probe receiver. Wind-driven surface waves were 
generated by a multi-fan blower placed at one end of the cavity. The 
relationship between the resulting fluctuation of the acoustic signal 
and the surface fluctuation above the receiver was investigated. The 
results are consistent with the simple theory at low surface-wave fre- 
quencies. Anomalous behavior was observed for the higher fr equencies 
of the surface wave spectrum. Possible mechanisms for this inconsistency 
are postulated and discussed. 
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